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Abstract. Dirac operators in non-trivial topology backgrounds in a finite box are re- 
viewed. We analyze how the formalism translates to the lattice, with special emphasis 
on uniform field backgrounds. 



I. Introduction 

Most of the numerical work performed within Lattice Gauge Theory makes 
use of periodic boundary conditions. This has the advantage of preserving 
translational invariance and homogeneity of the lattice. When the box size 
becomes large in physical units, correlation functions become independent 
of this size and boundary conditions. For smaller boxes and/or other observ- 
ables the effect of the boundary conditions is non-negligible. This fact has to 
be taken into account when comparing lattice results with the continuum. 

Topological aspects of the gauge fields are directly related to bound- 
ary conditions. There are well-known difficulties in extending topological 
notions to the lattice. For example, the space of three-dimensional (3D) 
non-abelian gauge transformations (or of 4D gauge fields) is connected on 
the lattice and disconnected on the continuum. In any case, it is more 
appropriate to relate topological properties of the lattice to those of contin- 
uum fields defined on the torus. For example, in trying to study individual 
instantons on the lattice there are, not only order a corrections, but also 
finite size corrections whose origin is in the continuum: there are no 4D 
self-dual solutions of topological charge 1 on the torus (without twist). 

The topological properties of gauge fields are directly connected with 
the spectrum of the Dirac operator through index theorems. Within Lattice 
Gauge Theories important progress in this respect has arisen lately from 
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considerations of Ginsparg- Wilson symmetry |0], domain- wall fermionsQ 
and the overlap method^]. It is therefore important to re-examine the 
study of the spectrum of lattice Dirac operators for non-trivial topology 
gauge field backgrounds. According to our previous reasoning we should 
consider the case of gauge fields in a box. In a recent paperQ the case 
of uniform (constant) field strength gauge field configurations in 2 and 4 
dimensions has been analyzed in detail (See Ref. || for an early work in 
this context). In this talk I will explain and develop certain aspects with 
greater depth. For simplicity and lack of space I will mostly concentrate in 
the two-dimensional case. We refer the reader to Ref. [|j for other aspects 
not covered in this paper and a more complete list of references. 

II. Two dimensions: Continuum 

Charged matter fields on a 2-torus are sections of a U(l) bundle. Without 
sacrificing rigor we can view them as complex functions ip(x) defined on 
the plane and periodic up to gauge transformations: 

ip(x + ei) = fli(x)ip(x) 

ip(x + e 2 ) = Sh(x)il>(x) . (1) 

where e\ = (Zi , 0) and e 2 = (0,fa) and U are the torus lengths. The U(l) 
fields fti(x) = exp{iu>i(x)} are the transition functions, and must satisfy 
the following consistency conditions: 

(u)i(x + e 2 ) — w\(x)) — (u}2(x + ei) — w 2 (^)) = 2irq . (2) 

where q is an integer characterizing the topology of the bundle. In math- 
ematical terms this is the first Chern number of the U(l) bundle. The 
physical interpretation of this integer can be deduced by considering abelian 
gauge potentials on this torus. They must satisfy: 

Aj(x + ei) = Aj(x) + djUJi{x) (3) 

Now the total flux of the magnetic field B = eijdiAj is given by 

J dx\dx2B = J dx2 c?2Wi(x) — / dx\ d\uj2(x) = 2irq (4) 

The explicit form of ooi{x) is a matter of gauge choice. A convenient widely 
used choice is: u>i(x) = nqeijXj/lj. A particular gauge field satisfying these 

boundary conditions is Af'\x) = —^eijXj, where F = 2^2 and A = l\h is 
the area of the torus. The corresponding magnetic field strength is constant 
and equal to F. On a torus the field strength is not the only gauge invariant 
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quantity, one has also the Polyakov lines winding around non-contractible 
loops. Hence, given a pair of real constants (defined modulo 2ir/li), 
one can construct a whole family of gauge inequivalent gauge potentials 
A^ v {x) = A^'{x) + Vi having the same constant field strength F. The 

fields play an important role in parameterizing the space of gauge 
fields compatible with the boundary conditions. Indeed, one can make use 
of Hodge theorem to write an arbitrary gauge field as: 

A = Af ] + Vi + di<t> + e^djh (5) 

where both 4>(x) and h(x) are real periodic functions on the torus. We see 
that Vi parametrises the harmonic 1-forms on the torus. Setting <f> = we 
have the expression of the gauge field in Coulomb gauge. 

Coming back to matter fields, let us denote by H. the space of func- 
tions defined on R 2 satisfying the boundary conditions Eq. [|. This has the 
structure of a Hilbert space. Therefore, given a complete set of commuting 
operators one obtains an explicit representation of this space. Covariant 
derivatives Di = d- L — %A{ with respect to compatible gauge fields are 
anti-hermitian operators defined on 7f, A natural basis is obtained by 
selecting the covariant derivatives with respect to the uniform field strength 
potentials: 

Bf = 9i - lA f > (6) 
They define a Heisenberg algebra: 

[Df\D^] = -rF (7) 

One can choose (to diagonalize) one the two operators to provide a natural 
representation of the space H. However, for \q\ ^ 1 this operator alone does 
not provide a complete set. There is an additional pair of operators 
commuting with Z?j and satisfying: 

K^R® =exp{^}K^K^ (8) 

defining a Heisenberg group. The boundary conditions imply (K^) q = I 
and the space decomposes into a direct sum of q subspaces. The presence 
of these q spaces corresponds to the finite degeneracy of Landau levels. In 
infinite space (q — > oo) the degeneracy becomes infinite and the relation 
leads to a new Heisenberg algebra relation. This is the quantum version of 
the canonical relations for the Hamiltonian system of a particle in a plane 
in an uniform magnetic field {P i = Pi± ^ € ij x j define the two commuting 
Poisson bracket algebras). 
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Our considerations lead to the following parametrization for an arbitrary 
element of Ti: 

ip(xi,x 2 ) = exp{nrq^-^-} ^ exp{2iri^(n + sq)} h n (x 2 + n + qS l 2 ) 

(9) 

where h n (y) are g arbitrary functions of a single real variable y. This yields 
a dimensionally reduced description of our Hilbert space. The price to pay 
is that, in general, differential operators acting on ip{x) map onto integro- 
differential operators acting on h n (y). However, covariant derivatives with 
respect to gauge potentials depending only on the variable x\ preserve their 
differential operator character. 

As a particular case one can consider the Dirac operator D in the back- 
ground of a uniform magnetic field. The spectrum consists on q-dimensional 
spaces of eigenvectors corresponding to the eigenvalues and ±i^/2Fp, 
where p is a positive integer. We see that the eigenvalues only depend on the 
area of the torus and not on its shape. The corresponding eigenfunctions 
in the h n representation are also universal, given in terms of harmonic 
oscillator eigenfunctions. Notice, however, that shape parameters enter, 
through Eq. |9], in expressing them in the standard basis. All these results 
hold for any value of u, since (for q ^ 0) covariant derivatives di — lA^ v (x) 
are unitarily equivalent. 

III. Two dimensions: Lattice 

Now let us examine the situation on a finite L\ x L 2 lattice. Transition 
functions do not appear in this formulation, so it seems at first unclear 
how one can make contact with the continuum formulation. However, if we 
consider compact U(l) lattice gauge fields in the standard way (Ui(n) = 
exp{— iAi}), it turns out that a unique decomposition similar to Eq. || holds 
on the lattice: 

Ai = +Vi + Af(f> + eijAjh mod 2nZ (10) 

where A 1 * 1 are the forward/backward lattice derivatives, and Vi are constants 
defined modulo 2ir/Li. The functions 4> and h are periodic on the lattice 
and uniquely defined up to an additive constant. Finally is the lattice 
version of the constant field strength fields: 

A (0) = -|Vi(n) (11) 

where Ii(n) = n\ for 1 < n 2 < L 2 and Ii(rii,L 2 ) = (L 2 + l)n\, and a 
similar definition for I 2 . The lattice constant field strength F = with 
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A = L1L2, contains information of the topological sector q. Apparently 
we have managed to uniquely define a topology of the lattice fields. The 
problem is that the decomposition Eq. [l(] is singular for fields having any 
plaquette on the lattice equal to —1. 

The lattice constant field is periodic. Replacing Ij by raj would lead 
to the naive non-periodic lattice discretization of the continuum field. The 
difference can be interpreted by considering that the lattice gauge potential 
includes the transition function Wi(n) = rrqeijUj /Lj. 

Lattice fermion fields ip{n) are elements of a complex vector space of 
dimension 2A. Usually one views these fields as being periodic on the lattice. 
However, fermion fields always appear in the action coupled to the gauge 
potential through covariant derivatives. Given the decomposition of the 
gauge potentials Eq. one can alternatively consider the non-periodic 
gauge fields Af 1 explained in the previous paragraph and non-vanishing 
transition functions. In this case, the fermion fields are required to satisfy: 

tp{n + e%) = exp{z7rg e^nj} ip(n) (12) 

where take the same form as in the continuum but expressed in lattice 
units (a=l). One might choose an appropriate basis on this space by requir- 
ing that a complete set of commuting operators (matrices) are diagonal. An 
important class of operators is given by the covariant shift operators %. All 
versions of lattice Dirac operators are elements of the algebra generated by 
these operators and the Dirac matrices. As in the continuum case we will 
attribute a special role to the covariant shift operators for uniform fields 
T^°\ They satisfy: 

T/ >7f>=r ? T 2 (°>7f> (13) 

where £ = exp{^} This relation defines a Heisenberg group. It has been 
studied in connection with twist-eaters §|]-@. One knows that the group 
acts irreducibly on a vector space of dimension A/(gcd(q,A)). This implies 
that for q = 1 one can fix a basis of the vector space by diagonalizing any 

one of the two 7^ . If gcd(q,A) > 1 there are operators commuting with 
these and allowing to define a complete set of commuting operators. If q 
divides Lj the formulas follow closely the continuum case. There is at least 
a q-fold degeneracy of all eigenvalues of any lattice Dirac operator. We refer 
to 01 for details. 



For simplicity, let us restrict to the q = 1 case. The condition Eq. |13 
has to be supplemented with the value of the Casimirs: 

(7; (0) )- 4 = exp{- i ^}I (14) 



All pairs of A x A matrices satisfying Eqs. p!3}-|l4] are equal up a similarity 
transformation. In our case, we have (pi = Avi, so that the spectrum only 
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depends on up to integer multiples of 2ir/A. This contrasts with the 
continuum where the spectrum is independent of Vi. 
In the basis which diagonalizes we have: 

T} 0) = Q = exp{-«h./A} diag{l, £, £ 2 , . . . ^} (15) 
_ p w ith Pjk = 5kj+i exp{— vfa/A} (16) 

Thus we see that, irrespectively of the lattice Dirac operator we are using, 
the spectrum does not depend on the shape of the box but only on the area 
A (and Vi mod 2ir/A). Indeed, the spectrum also coincides with that of a 
twisted reduced model for the group J] (A) §§ at weak coupling. 

Now we will investigate the properties of the spectrum of the Wilson- 
Dirac hamiltonian HwD(M,r) = ^^,DwD{M,r) where Dwu(M,r) = M + 
Dn + rW. The naive lattice Dirac operator Dn and Wilson term W are 
given by: 

i 

W = \^- + (18) 

i 

We explicitly indicate the dependence of the Wilson-Dirac operator on the 
mass M and the Wilson parameter r. 

In analyzing the spectrum of H\y£>(M,r) it is useful to make use of 



symmetries. They must correspond to symmetries of the algebra Eq. 13 
For example the transformation 

7^>— >7^>; T^^T^ (19) 



preserves Eq. |13j. However, it can only be an exact symmetry if it preserves 
the Casimirs as well, and this only occurs for fa = 0, n. In this case, it can 
be combined with a similar rotation in spin space (7, — ► e^Tj) to produce 
a unitary matrix U commuting with H\v£>(M,r). The operation U can be 
interpreted as a ir/2 rotation of space, which happens to be symmetry even 
when L\ ^ L2. Indeed, it generates a finite group of 4 elements U 4 = I. 

For even A one can construct additional symmetry operations associated 
to 7^ — > —T^°\ Using them one can show that Hw£>(M,r), Hwd(M + 
4r, — r) and —H\yd{—M — 4r, r) are unitarily equivalent. 

In Ref. p| formulas were given that allow the computation of the eigen- 
values of Hw£>(M,r) to machine precision. One particularly interesting 
aspect is the balance between the number of positive (N+) and negative 
(iV_) eigenvalues of Hwd(M, r). Obviously we have N + + iV_ = 2 A, while 
X = N- — A provides a lattice definition of the index. It is precisely the 
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index of the Neuberger (overlap) operator. For large positive or negative 
values of the mass M it is easy to show that N+ = and the index 
vanishes. As we decrease the mass towards negative values some eigen- 
values might move from positive to negative or vice versa generating a 
non-zero index. Jumps of the index take place at values of the mass M 
for which there exist ij)^ satisfying Hwd(M, r)ipj^ = 0. Multiplying by 
73 the previous equation, one sees that m ust be an eigenvector of real 

eigenvalue (equal to —M) of Dwd(0, r). Furthermore, the sign of V'jgTsV'jw 
determines whether the index increases or decreases at this point. Thus, 
to analyze the index of H\vd{M^t) it is enough to determine the real 
eigenvalues of Dwe>(0, r) and the corresponding eigenvectors. For r = 1 one 
can apply similar techniques to those used in Ref. § to derive expressions 
that allow the computation of this real eigenvalues to machine precision. 
For A > 3 there are always four real eigenvalues, two in the interval 
[—2,0) and two in the interval (—4,-2] (For even A there is symmetry 
around M = —2). As we decrease the value of M from positive values 
the index follows the sequence 0— > 1 ^ — > — 1— > 0. The physical 
region at which the index takes its continuum value (equal to one) has an 
upper edge given by -0.63397,-0.58578,-0.29449,-0.12210,-0.0311755 
for A = 3,4,10,25,100. For large torus sizes it approaches —Fr/2. The 
lower edge is —0.94203 for A = 3 and decreases very fast towards —2 as 
the area increases. One can monitor the M dependence of the smallest (in 
absolute value) eigenvalue of Hwd{M, r). The interval (—4, 0) is split into 4 
regions, one for each of the eigenvectors of real eigenvalue tp^. Within each 
region the behaviour is essentially linear (in M) with a slope determined 

by ^ Si i^m- 

Other aspects are studied in Ref. [|[]. For example, we numerically com- 
puted the eigenvalues of Neuberger's operator in a uniform background 
gauge field, for several torus sizes up to 24 x 24 and several values of M. 
Eigenvalues and eigenvectors are quite close to the continuum results. 

IV. Four dimensions 

Much of the previous construction generalizes to the case of non-abelian 
gauge groups in 4 dimensions. For SU(N) gauge groups the topology of 
gauge fields is indexed by the instanton number Q (and twist sectors). To 
make connection with the two-dimensional U(l) formulation of the previous 
sections, we make use of the fact that (except for SU(2) and odd Q) it is 
possible to construct transition functions living in the maximal abelian 
subgroup U(l) JV " 1 Q- 10| ■ Then the boundary conditions for the fermion 
fields apply for each of the color components independently. This reduces 
the problem from SU(N) to a U(l) field in 4 dimensions. The fluxes of these 
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U(l) fields over the 2-planes encodes the information of the topology of the 
original field. Choosing appropriate transition functions and a basis of the 
4 dimensional space, the problem decouples into two 2-planes. This allows 
to obtain the spectrum of the continuum Dirac operator and the lattice 
naive Dirac operator in terms of the two-dimensional result. The Wilson- 
Dirac and Neuberger operators, however, couple the 4 dimensions in a more 
complicated way. We refer to Ref. |3j for details. 
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